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Abstract 

The heat transfer model for a one-dimensional supercooled melt during the 
final stage of solidification is considered. The Stefan problem for the de- 
termination of the temperature distribution is solved under the condition 
that (i) the interface approaches the specimen surface with a constant veloc- 
ity V; (ii) the latent heat of solidification linearly depends on the interface 
temperature; (iii) all the physical quantities given at the phase boundary are 
presented by linear combinations of the exponential functions of the interface 
position. First we find the solution of the corresponding hyperbolic Stefan 
problem within the framework of which the heat transfer is described by the 
telegraph equation. The solution of the initial parabolic Stefan problem is 
then found as a result of the limiting transition V/Vh — > (Vh — » oo), where 
Vh is the velocity of the propagation of the heat disturbances, in which the 
hyperbolic heat model teds to the parabolic one. 

Keywords: Solidification, Final transient, Heat transfer, Stefan problem, 
Telegraph equation 

1. Introduction 

The process of rapid solidification is a well established method for the 
production of metastable solid states of different nature making it possible to 
study new mechanisms of crystal growth and produce materials with radically 
new physical properties 
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In rapid solidification experiments very high velocities of phase inter- 
face can be reached. Such conditions occur during solidification of the un- 
dercooked melts or recrystallization after pulsed-laser irradiation of a solid 
surface. When the interface velocity reaches some critical value diffusion- 
temperature field in the bulk of both phases can significantly deviate from 
local equilibrium 0, [sf] . In this case both the diffusion and heat fluxes are no 
longer defined by the classical Fick's and Fourier's laws relating the diffusion 
and heat fluxes correspondingly to the gradients of a solute concentration 
and temperature. The simplest generation of Fick's and Fourier's laws tak- 
ing into account the relaxation to local equilibrium in the diffusion and the 
heat field is given by the Maxwell-Cattaneo model and leads to the hyperbolic 
transport equations jsj. 

In the past two decades a great body of studies, devoted to the local 
nonequilibrium heat and mass transport during rapid solidification, has been 
executed @, S 3, 0, @, 0, S 0, (U (3 01 Q IS (3 E3 [3 • The numerical 
estimates show that under conditions of experimentally achievable interface 
velocities local equilibrium is only disturbed in the diffusion field, while the 
heat field can be described in the local-equilibrium approximation within the 
scope of the conventional parabolic heat conduction model 0, 0] 

The currently existing analytical models of the directional solidification 
processes usually consider the initial transient and the motion of the planar 



front far from the boundaries of a system [19|, [20j, [21]. The investigation of 
the final transient is practically absent. Meanwhile, besides a purely aca- 
demic problem there exists considerable practical interest as well because 
the final study of solidification process influences the formation of the sur- 
face layer of the materials, their surfa ce phys ical-chemical characteristics and 



the distributions of different defects 22, 23, 24 



The final transient of a binary melt solidification has been analytically 
considered by Smith et. al. in the work 25) (also see 26), p. 278). In the 



local equilibrium approximation the authors have calculated the terminal 
solute distribution of a formed solid. However the found distribution has 
divergence at the surface of the specimen that may be caused, among other 
reasons, by neglect of the temperature changes during the interface motion 
near the surface. Because of the thermodynamic relationship between the 
interface temperature and the solute concentration, resulting from the phase 
diagram, such changes can lead to further solute redistribution in the bulk 
of the specimen. In this connection it is of interest to initially investigate 
the evolution of the temperature field of the solidifying pure melt during the 
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final transient. 

In general, the problem of the determination of the temperature field dur- 
ing solidification of a melt is known as the Stefan problem. It consists in solv- 
ing of the heat conduction equation for a temperature T of each phase with 
the boundary conditions at the moving interface and in the one-dimensional 
form is written as 

dT L . s d 2 T L , s 

T L = T S , (2) 
q L -qs = QV, (3) 

where indexes L and S are respectively related to the liquid and the solid 
phase, p is the density, c p is the specific heat and A is the thermal conductivity 
(for simplicity, all material characteristics are assumed to be constant and 
identical within the phases and at the interface). Equation (j2j) represents 
the condition of the continuity of the temperature across the interface and 
equation (|3]) defines the condition of a heat balance at the interface, where 
q^s — —^dT Lt s/dx is the heat flux, Q is the latent heat of solidification and 
V is the velocity of the interface. In addition, at the surface of the specimen 
the boundary conditions must be given. 

From the thermodynamic point of view the velocity V is determined by 
the undercooling of the interface AT = T m — T, where T m is the equilibrium 
temperature of solidification and T is the interface temperature, V = /(AT), 
so that at AT = 0, V — /(0) = 0. For so-called normal crystal growth it is 
assumed that 

V = a«AT (4) 
where \i is the kinetic coefficient characterizing atomic attachment kinetics 



at the interface [19j. When the undercooling of the interface AT is large 
enough, the relationship (j4]) can no longer hold. This takes place when the 
melt is initially supercooled. In a number of experiments with pure metals, as 



well as by means of molecular-dynamic simulation [27], 128|, [29|, [30[ it has been 
shown that at the beginning the growth velocity increases with increasing 
undercooling AT reaching a maximum value and then above some critical 
undercooling AT*, V is practically kept constant in some region AT (also see 



19| . p. 18). In what follows we shall presume that the function V = /(AT) 



has such properties. 
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When the undercooling is large enough the amount of latent heat released 
on the solidification front may prove to be deficient to heat the interface to 
the temperature T m . In this case the undercooling at the interface will be 
different from zero during the whole solidification process, AT ^ 0, and 
the solid phase will reach the surface of the specimen with a finite velocity, 
Vf = f(ATf) 0. When the interface moves in the near surface region 
the undercooling AT, generally speaking, will change, at the same time also 
changing the growth velocity V. However, assuming that Vf is high enough, 
due to high undercooling and AT changes within the region for which 
velocity V depends slightly on AT, one can consider that the interface moves 
near the surface with approximately constant velocity V equal to Vf. It 
should be noted that there exists a substantial distinction from the situation 
in the initial transient when the interface velocity changes from zero to a 
steady state value. 

The behavior of the temperature T, the heat fluxes qi,s and the latent 
heat Q at the moving interface essentially affects the evolution of the thermal 
field in the bulk of the phases. Giving different models of their behavior at 
the interface, one can consider various models of the solidification process. 
In the present study we consider the exactly solvable model of solidification 
within the scope of which any physical quantity F (the temperature, the heat 
flux, etc.) is given at the interface as 

F(x,t)\ t =m = A { F) + Af\~^ + A<Fe-T*'* • • • , (5) 

where t = t(x) determines the path of the interface position and the coef- 
ficients A n and the powers of the exponents 7 n must be defined from the 
phase boundary conditions and the boundary condition at the surface of the 
system. It is worth noting that an expression of similarly type has been 
obtained for the solute distribution in the work 25 . 



As regards the latent heat of the solidification it is normally assumed to 
be constant and is empirically defined for the equilibrium temperature T m 
as Q = Q m = kT m (k = const) ((HI, p. 185). If the undercooling of the 
interface is high enough, it is reasonable to define the variable latent heat as 
Q = kTi = Q m — kAT that will be used in what follows. 

Thus for the determination of the temperature field within the scope of 
the given model we seek the solution of the one-dimensional heat conduction 
equation (JTJ with the boundary conditions (j2J), © and (jSJ) when the interface 
approaches to the surface with a constant velocity V and the latent heat 
linearly depends on the undercooling AT. 
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As has been noted above, the temperature field in a rapid solidifying melt 
can be considered within the scope of the parabolic model ([I])-©- However to 
obtain the solution of interest we shall initially consider an auxiliary problem, 
namely, the corresponding Stefan problem for the hyperbolic heat conduction 
equation (see, for example, j5|,[l2[). It turns out that the boundary conditions 
given at the moving interface are more simply taken into account within the 
scope of the hyperbolic model. As is known, the hyperbolic model of heat 
conduction, based on the telegraph equation, gives the finite velocity of the 
propagation of the heat disturbances in matter Vh and is reduced to the 
parabolic model ([[]) in the limit Vh — > oo [8]. The idea of the work is that 
initially the solution of the hyperbolic Stefan problem with arbitrary ratio of 
the velocities V/Vh < 1 is solved and then the limiting transient, V/Vh —> 0, 
to the solution of the parabolic problem is executed. 

The work is organized as follows. In Sec. 2 the hyperbolic Stefan problem 
corresponding to the boundary problem flTJ)-(j3]), (jSJ) is considered. The solu- 
tion of the telegraph equation is found by the Riemann method within the 
scope of which the boundary conditions given at arbitrary moving bound- 
ary are automatically taken into account. On the basis of this solution the 
heat fluxes and the temperature fields both in the liquid phase and in the 
near interface region of the solid are determined. The subsequent limiting 
transition V/Vh — ► and the solution of the parabolic problem are given in 
Sec. 3. The conclusion is presented in Sec. 4. The Riemann method and its 
application to the presented problem are contained in the Appendices. 



2. Hyperbolic model 

The hyperbolic model of the heat conduction starts from the Maxwell- 
Cattaneo relaxation equation for the heat flux [§]. In the one-dimensional 
form this equation for the liquid phase is 

qL + T ^T = ' (6) 

where r = a/V H is the time of the relaxation of the heat flux to its local 
equilibrium value defined by the Fourier's law and a = \/pc p is the thermal 
diffusivity. 

Equation (jSj) in combination with the conservation law 
9Tl dq L 

PCp — (7) 
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gives rise to the hyperbolic transport equations 
d 2 T L dT L d 2 T L 

d 2 q L , dq L d 2 q L 
T d¥ + dt ~ a dx 2 ■ [ ' 

The equation of the type (jSJ) or Qj is known as the telegraph equation. At 
r — > (or Vh — > oo) the equations (jSJ) and fl5]) are reduced to the parabolic 
heat conduction equation (jTJ). 

Now let us consider the supercooled pure melt initially occupying the 
half-space x > 0. The planar front of solidification forms in the infinitely 
removed region at t = — oo and moves in parallel to the specimen surface fixed 
at x = 0. As it has been noted in the Introduction when the undercooling 
is large enough the interface will move in the near surface region with the 
approximately constant velocity V — Vf along the path x + Vt = 0. In this 
case the region occupied by the melt in the final stage of the solidification 
process is given by the inequality ^ x ^ — Vt (t ^ 0). Therefore in the 
plane (x, t) the liquid phase occupies the region x + Vt ^ 0, x ^ 0, t ^ 0. At 
the interface the condition of the heat balance Q holds 

(Ql - qs)\x+vt=o = -VQ\ x+Vt =o , (10) 

Now we consider the heat flux qi in more detail. Introducing dimension- 
less variables t/r, x/tVh in the equation (Q, one obtains 

dPq L dq L = d 2 q L 

dt 2 dt dx 2 ' 1 ] 

where the former notations (x, t) have been used for new variables, q~L = 
QlUQtuVh) is a dimensionless heat flux. The boundary condition ( TTOj) in the 
dimensionless form is written as 

(Ql ~ qs)\x+cd=o = -aQ\x+at=o , (12) 

where q$ = qs/iQmVn), Q = Q/Q m and a = V/Vh is the dimensionless 
parameter. In addition, we assume that at the surface the equality should 
be fulfilled 

q L {xt)\ x=0 = (*<0), (13) 

expressing the condition of the absence of the heat flux through the surface. 
Finally, the solution of ffTT]) is sought in the near surface region at X = 



6 



x + at ^ 0, x ^ 0, t ^ occupied by the liquid phase while the solid 
occupies the region X ^ (see fig. lA.lb in Appendix). 

Now we consider the case of a < 1. Suppose that at the moving interface 
residing in an arbitrary point x near the surface at the moment t = —x/a 
the flux c[l and its time derivative dqi/dt are known 



U{xt)\ 



t=-x/a — Qoyx) 



dqdxt) 
dt 



\t=-x/a — Ql[%) j 



(14) 



where the functions qo(x) and qi(x) will be specified further. 

If the functions qo(x) and qi(x) are known the solution of the equation 
pip satisfying the conditions (Till) in the region X ^ at a < 1 can be found 
by the Riemann method [32| (for details see Appendix A) and has the form 



q L (xt) 
1 



x + t 

2 in- a r^ 



exp 



X 



12(1 - a)l 



x — t 
+ — )exp 

' 1 + a 



X 



2(1 + a) 



a.(x — t) 
1+a 



dxMxi)e- Xl/2a J (I v 7 ^ - ^i) 2 " (t + xi/a)'- 



X 
4a 



a(x+t) 
1-a 

a(x — t) 

1 + Q 



-t/2 



dxMx^e-* 1 ^ 



4 (§ ^{x-xrf-it + xr/af) 



a(x + t) 
1-a 



\/(x - xi) 2 - (t + xx/a)* 



(15) 



where 



a- 



g x (z) 



(16) 



and Jq(x) is the Bessel function of zero order. 

In accordance with what was said in the Introduction all the quantities 
given at the phase interface are represented by linear combinations of the 
exponential functions (jSJ). In particular, let <p(x) and if)(x) be given by the 
expansions 



<p(x) 
ip(x) 



A + A ie ^ ix/2 + A 2 e~^ x/2 ■ 
B + B ie ^ ix/2 + B 2 e~^ x/2 ■ 



(17) 
(18) 
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where constants 7„ ^ 0, A n and B n will be specified in what follows. After 
the substitution of (ITTj) and (1T8T) in (fi3|) and the calculation of the integrals 
(details in Appendix B), we obtain 



q( x t) = ^V^j^exp 



ln +) X 
2(1 -a 2 



+ 4+) 



cxp 



2(1 -a 2 



where the following notations have been introduced 
7n + a ± a/ a 2 7 2 + 2«7„ + a 2 ^ ; 



7 (±) 

in 



± _ An , „ #n . 

— n — ' 

2 V n 



a 



1 + «7„ 



/A, 



1-^(1"« 2 )- 
or 



(19) 

(20) 
(21) 

(22) 



Let us determine the parameters j n , A n and B n in such a way as to satisfy 
the balance condition ffl2|) and the boundary condition at the sample surface 



2.1. The determination of the parameters 

Now consider the boundary condition (fT3"j) . Taking into account that 
7o = 0, 5q = a, uq = a, = a ± a and using the equation ( TT91 . we have 
for arbitrary small t < 



x=0 



2aH , A (+) 



= exp 
+ aS _) exp 
+ A { 2 ] exp 



2(1 


-a 2 ) 


Tr" 


f W 


2(1 


-a 2 ) 


7* 




2(1 


-a 2 ) 



+ 



A[ +) exp 



7i W 
2(1 - a 2 ) 



+ 



a(+) ^2 ^ 

^2 ex P 2T + 

2 1 — ar 



0.(23) 



If all the powers of the exponentials are different then q(0, t) = can be 
only at A n = B n = 0. However if each exponential function appears in the 
equation (|23|) at least twice then this can lead to nonzero A n and B n . Bearing 
in mind this circumstance we determine 7„ so that the following equalities 
hold 



7 H = 7 M 

in in—1 



n — 1, 2, 3, . . . 



(24) 
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Table 1: The parameters of the equation (|T 



n 





1 


2 


3 


4 


In 
In 
In 




2a 



4a 


4a(3+a^) 


80(3+0^ )(l+a 2 ) 


8a(l+a^)(a 4 +10a 2 +5) 


1-a 2 

8a 
1-a 2 

2a 


(1-a 2 ) 2 
2a(3+a 2 ) 2 
(1-a 2 ) 2 
8a 
1-a 2 


(l-a 2 )3 
32a(l+a 2 ) 2 

(l-a 2 )3 
2a(3+a 2 ) 2 

(1-a 2 ) 2 


(1-a 2 ) 4 
2a(a 4 +10a 2 +5) 2 
(1-a 2 ) 4 

32a(l+a 2 ) 2 
(l_ a 2)3 



in which 7^1_\ (and respectively 7„_i) are considered to be known 0. Taking 
into account the notation (120]) and resolving the equation (124"]) in relation to 
7„, one obtains 



(1 - a 2 )( 7n ) 12 = 7 « ± A /a 7 a[2(l-a2) + a 7 ^ ) 1 ] . (25) 



At n = 1 and 7q + ' ) = 2a the equation (125]) gives 
4a 

7i = 



a 2 



The second value 71 = is the extraneous root of the equation (124]) at n = 1. 
After the determination of 71 the values 7^ appearing in ( [19]) can be found 
from the equation ( 120]) . Along a similar line one can obtain the values 7„, 
7n for n > 1. In Table H] these values are given for n < 4. As is seen from 
the table 7n ~ (1 - a 2 )"", 7 i +) ~ (1 - a 2 )" n 7 i~ } ~ (1 - a 2 )~ n+1 . The case 
of an arbitrary n is easily proved by induction using f ]25]) . 
Under condition f ]24"]) . the equation (]23l holds, if 

4 +) = ^0/2 + ^0 = 

4~ } = A /2-5 = A (26) 

Finally taking into account the equalities (126]) . the expression ( 1T9]) can be 
rewritten in the form 

q(xt) = ^M^ n+ix/2 ~ e" 7W2 ) exp [ J • (27) 



^he equation 7n = 7.^1 either has no the solutions or does not give the new ones. 
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2.2. The temperature field 

The temperature field in the liquid phase can be found in the same way as 
the heat flux has been defined. The resulting expression for the dimensionless 
temperature Tl takes the form 

f L (xt) = 4 +) + J2i a n )e ~ 7nX/2 + 4V 7 " +W2 } exp Lffoj > ( 28 ) 

where T L = pc p (T L — T m )/Q m . The constants aft^ 1 can be expressed in terms 
of the parameters determining T and dT/dt at the interface by the equations 
of the type (jlgjl-fllgjl and j2Q. 

Substituting the expressions for the flux ff27|) and the temperature fl28l) 
into the energy consideration law (jTJ) and equating the coefficients at the 
linear independent functions, one can express the constants On in terms of 
appearing in (127|) . The corresponding expressions will be given for the 
case of the parabolic model. 

2.3. The solid phase 

The heat flux q$ and the temperature Tg in the solid satisfy the equations 

d 2 q s dq s _ d*q s d 2 f s df s _ d 2 f s 

dt 2 dt dx 2 ' dt 2 dt dx 2 [ ' 

where q s = qs/(Q m V H ), f s = pc p (T s - T m )/Q m . 

For the complete determination of the temperature field in the liquid the 
interface boundary conditions ([2]) and f|T21) depending on the solid tempera- 
ture and the heat flux must be used. For their determination it will suffice to 
consider the solutions of the equations (|29|) in the region near the interface 
defined by the inequalities X>0, x + t<0 (see fig. lA.lb ). The solutions 
of the equations (1291) in this region can be obtained in the same way as for 
the liquid phase. The application of the Riemann method in the indicated 
region gives for the heat flux 

t) = £ e--/'^<-> exp +i<+» exp } , (30) 



q s {xt 



where 7^ are given by the equality fTSDl and the constants can be 
expressed in terms of the parameters determining the flux q$ an d its time 
derivative at the interface by the equalities of the type (I2T|) . 
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The expression for the temperature T$ is analogously written down as 



f s (xt) = £ e-W2{ a H exp + 5 ( + ) exp [_%_^_] }. (31) 



2(1 - a 2 )l n ^L2(l-« 2 ) 



3. The parabolic model 



As has been indicated above the transition to the parabolic model is 
executed by the limit a = V/Vh — > (V# — > 00). Using table [1] and the 
equation ( 125]) it is easy to show by the induction for any n that for small a 

7 n ~ 2n(n + l)a a -» (a^O) 

7^ _) ~ 2n 2 a (32) 
7 (+) ~ 2(n + l) 2 a. 

When the relationships fl32|) are fulfilled the expressions for the temperature 
Tl and the flux in the liquid phase are written down in the form 

f L (xt) = 4 +) + ]T r- Jn+1X /2 + e -^/2 }e ^[X/2 (33) 

q L ( xt ) = An+i{ e ~ ln+lX ' 2 - e-^ x ' 2 }e^ x/2 (34) 

The constants On in the equation (|28p for have been defined in such a 
way as to satisfy the conservation law ([7]) (see the end of section 2.2). 

From expression fl3"3"|) it is seen that the disturbances of the temperature 
field ahead of the solidification front propagate only over distances in the 
order of I < 2tVh/'^\~" > = a/V (in the dimensional variables). Therefore, 
if the interface is removed from the surface at the distance I ~ a/V, the 
surface still remains at the initial temperature T H. It is supposed, of course, 
that the constant interface velocity approximation holds over distances in the 
order of a/V from the surface. At x = and V\t\ ~ a/V in the expression 
(I3"3"j) one can neglect by sum (\X\ = \Vt\ ~ a/V) and write down for the 
temperature at the surface 

J-L\x=0 ZA ~ Oq , 



2 For example, for Ni a = 12 • 10 6 m 2 /s V ~ 20m/ s, a/V ~ 0.6^m. 
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where A = pc p (T Q — T m )/Q m < is the initial undercooling of the melt. 

Now let us consider the temperature field in the solid phase. The equal- 
ities f l30|) and fl3Tl) hold in the region between the straight lines x + at = 
and x + 1 = (see fig. lA.lb ). or in the dimensional variables, between the 
straight lines x + Vt — a; + Vfft = 0. At Vff — >■ oo the second line goes to 
the straight line t = 0,0<a;<oo and the region of interest to us will be 
given by the inequality — Vt < x < oo, spreading over the whole solid phase. 

The variable part of the expressions (130]) and ( I3T]) in the dimensional 
coordinates (x, t) is determined by the exponents 

e (n+l)V[x+(n+l)Vt]/a ^ e -nV[x-nVt]/a (V/V H < 1) . 

It is easy to see that at small t the terms containing the first exponent (are 
proportional to or 5rT ) with increasing x will indefinitely increase. 
In order to avoid such nonphysical behavior we put = = and 
introduce the notations A^ = A^. Turning back to the dimensionless 
coordinates, let us write down the expressions ( 15D|) and fT5Xj) at the small a 
in the form 



f s (xt) = a { s) + ^e-^ 2 e^ >x / 2 (35) 



(S) 



an 



U*t)=J2 Ai n S) ^ nX/2 ^ )X/2 > (36) 



where the constants ai + ^ have been determined so that the conservation law 
((7j) is obeyed. It is easy to check that the expressions fl33|) - fl36|) satisfy the 
heat conduction equations 

dT LS _ d 2 T LS dg LS _ d 2 q LS 

dt dx 2 ' dt dx 2 

and the Fourier's law is fulfilled, qi,s = —dTL t s/dx. 

3.1. The temperature field 

For the determination parameters appearing in equations (l33|) - (!36|l we 
use the condition of continuity of the temperature across the interface §Z§. 
The detailed calculations are given in Appendix C. The final expression for 
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1 2 3 4 0,0 0,5 1,0 1,5 

xV/a xV/a 



Figure 1: The temperature profiles for different moments of time. The dashed lines corre- 
spond to the temperature distribution in the solid phase, the solid lines correspond to the 
liquid phase; t = a/V 2 . 



the temperature of the liquid phase can be represented in the dimensional 
coordinates (x, t) as 

;i + 6A) 



f L (x, t) = A + 1 _^ J2 C n {e nV ^ +nVt ^ a + e-^^-nVtj/a j (37) 

1 

{2k - 1 - 6) 



C^U C n = b^H (9h _^_ h) (n>2); 

k=2 

(0 < x < -Vt, t < 0) : 



where 6 = Tq/T m and Tq = Q m /pc p . For metals the dimensionless parameter 
6 varies through the range < b < 1. For example, for Ni, T m = Tq = 



<2m/pc p = 397K and 6 = T Q /T m = 0, 23 [id 



Similarly one can write down for the solid phase 

M ' ; 1-6 1-6 ^ 2n + 1-6 ' K J 

It is easily seen that each term in the brace (J3"7]) represents the superposi- 
tion of two heat waves propagating in the mutually opposing directions with 
the velocity nV. 

3.2. Numerical results 

Figures [T]l2] present the temperature profiles obtained from the equations 
fl3ZD, d3SD for A = -2.5 (supercooling) and 6 = 0, 23. 
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b=0.3 



b=0.25 



b=0.2 



-2,2 







2 
xV/a 



3 



4 



Figure 2: The temperature at the interface depending on the distance to the free surface 
for different values of the parameter b = Tq /T m . 

The temperature curves for some moments of time are shown in figure 
CD The dashed lines are the temperature distributions in the solid phase, the 
solid lines give the temperature field in the liquid phase. As is seen from figure 
[T^i when the interface is relatively far from the surface (t/t = —3, t = a/V 2 ) 
the temperature of the solid phase is constant, and the temperature of the 
liquid phase falls to approximately the initial temperature of the melt T 
(T L \ X=0 « A = (T — T m )/Tq = —2.5) at the surface. When the interface 
moves close enough to thermal isolated surface (figure [lb), the released latent 
heat gives rise to the gradual heating of both the liquid phase and the near- 
interface region of the solid. 

In the figure [2] the dependence of the interface temperature on the param- 
eter b is shown. From the figure it is also seen that the interface appears on 
the surface (at x = 0) in the supercooled state, |T — T m \/TQ > 1, providing 
high final velocity of the solidification processes V — Vf. 

4. Conclusion 

In the given work we have considered a one-dimensional model of the heat 
conduction in the supercooled melt during the final transient. Three main 
assumptions underlie the model. Firstly, it is supposed that the interface 
approaches the system surface with a constant velocity. Some reasons for this 
assumption are provided by a number of experimental results and molecular- 
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dynamic simulation 27], |28|, [29|, |30|, [l9| showing that when the undercooling 



of the melt is large enough the interface velocity can slightly depend on 
the temperature. The second supposition assumes that the latent heat of 
solidification linearly depends on the interface temperature. Finally, it is 
supposed that the physical quantities of interest (the temperature, the heat 
flux, etc.) given at the interface are presented by linear combination of the 
exponential functions of the form (J5j), the parameters of which are determined 
as part of the general solution of the problem. 

Within the scope of the model the exact solution of the one-dimensional 
Stefan problem ([I])-®, ([3} defining thermal distribution in the system when 
the interface moves near the surface has been found. To this end, initially, 
the corresponding hyperbolic Stefan problem has been considered within the 
framework of which the heat transfer is described on the basis of the telegraph 
equation. The telegraph equation for the heat flux and the temperature in 
both the liquid phase and the near interface region of the solid has been 
resolved by the Riemann method. Further we have used the fact that in the 
limit a = V/Vh — > the hyperbolic heat model is reduced to the parabolic 
one. Taking into account this circumstance and executing the limiting tran- 
sition a — > in the expressions for the fluxes (127|) . (1301) and the temperature 
(|28|) . (I3T]) the thermal distribution in the sample during the final stage of 
solidification has been obtained. 

In conclusion it should be noted, that the given approach allow us to 
consider also other models of the solidification process differing from model 
(5). It is likely that the interface boundary conditions in the form of the 
superposition of exponential functions are the only ones for which the exact 
solution exists. On the other hand the solution for the heat flux (15) is 
written down for arbitrary boundary conditions (arbitrary if and iji) and 
opens up the possibility of numerical simulation. 
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Figure A.l: (a)The figure to the Riemann method, (b) The solution of the equation ([TTj) is 
sought in the region of the liquid phase x + at < 0. (c) The solutions of the equations ([29| 
are sought in the region of the solid phase situated between straight lines X = x + at = 
and x + t = . 



Appendix A. The Riemann method 

Let it be required to find the solution of the linear hyperbolic equation 



d 2 q dq d 2 q 
dt 2 dt dx 2 



(A.l) 



satisfying the initial conditions given at the curve T : t = g(x) (see figure 

E3k) 



q\t=g(x) 

dq 

dt t=g(x) 



q x {x) . 



The substitution q = e t / 2 u makes it possible to lead equation flA.lj) to a 
more simple form 

d 2 u d 2 u 1 



8x 2 dt 2 + a u ° 

with the initial conditions 



u \t=g(x) 

du 



dt 



t=g(x) 



q (x)e^ 2 = ^(x) 
{\qo + qi)e 9{x)l2 =^{x). 



(A.2) 

(A.3) 
(A.4) 
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The characteristics of equation (1A.2D are the straight lines x ± t = const. 
According to the Riemann method [32J if the characteristics go through the 
point M and intersect with the curve T at the points P and Q , then the 
solution of equation flA.2j) at the point M can be represented as 



u(M) = -iu P + u Q )- 



1 f ( du , du \ / dv . 9« , \ . . 

2 7 u fe^ 1 + df^ 1 ) - u KbV^ 1 + dT^ 1 ) (A ' 5) 



PQ 



The integral in (1A.5|) is taken along the curve V from P up Q and «p and 
uq are the values of u, taken at the points P and Q. The Riemann function 
v(M, Mi) for equation (IA.3p has the form 



v(M, Mi) = Jo (| v 7 ^ - xO 2 - (t - ^ 



(A.6) 



where Jo(x) is the Bessel function of zero order and du/dx is calculated along 
the curve as 



du 
dx 



t=g(x) 



<p[(x) - xpx(x)g'(x) 



(A.7) 



The Riemann method for arbitrary linear hyperbolic equations can be found, 
for example, in (32[. 

Now consider the solution of equation ( 1A.3I) in the region x ^ 0, t ^ 0, 
X = x + at < 0, when the initial data are given at the straight line t = —x/a 
(see Figure lA~Tb ). Instead of ( 1A.3j) and (1A.4j) we have 



u 
du 
~dt 



t=—x/a 



-x/a 



q {x)e 
1 1 



-x/2a 



Qo + qi)e 



-x/2a 



(A.8) 
(A.9) 



If the point M has coordinates (x, t), so it is easy to show that the points P 
and Q have the abscissas respectively equal to 



Xp 



a(x + i) 
1 — a 



a(x — t) 
1 + a 



(A.10) 
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Consider the integral term in equation (1A.5[) . Using equations f1A.7[) - (1A.10n 
and the fact that along the pathway of integration dt\ = —dxi/a, one has 



a.(x — t) 

1 + Q 

1 ' dx.e-^iv^+^xM-^-^-)) , (A.ll) 

l. \aOXi Ol\J > ti=-xi/a 



1—a 



where the notations are introduced 
ip{x) = q (x) , 



111 oP 1 

^0) = o?o(^) q'o(x) 2~9l(^) 

2 a a 



Furthermore using the Riemann function ( 1A.6I) . it can show that 



a dxi dti ) 



1 dv dv \ v -\M-r -./•)'- ir ; ./• /n ' 2 



h^-xi/a 2a y/(x - Xi) 2 - (t + x~Jaf 



.(A. 12) 



Finally, after substitution of integral (lA.llj) into equation (1A.5j) and using 
the equality q = e~'/ 2 w, one obtains the solution of the starting equation 
(lA.ip . with added conditions ( !T4l) . in the form represented by ( fl5l) . 

Appendix B. The calculation of the integrals 

Substituting equations ( ITT]) and f|T8|) into f|T5|) we have 

q L (x,t) = ^J n (x,t), (B.l) 

where 



x,t) 

T 



+ 



-B n jj ll ^ -\- A n J^ ^ -\- 
~aj n (x + t) + X 



cxp 



2(1 



a) 



+ exp 



a^ n (x — t) + X' 
2(1 + a) . 



(B-2) 
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and 



i e -t/2 



a(:E — t) 

1+a 



eteie-^/^Jo^v^-Zi) 2 - (t + Xi/a) 2 ) ; 



4 2) = ^e"'/ 2 
4a 



*(j;-t) 
1+a 



d Xl e- Xl/25 > 



(B.3) 

J£(|V(x-Xi) 2 -(t + Xi/a) 2 ) 



1-a 



yf{x — xi) 2 - (t + xi/af 



a 



1 + «7 r , 



(B.4) 
(B.5) 



T/ie calculation of J, 



(i) 



Making the substitution in the integral (1B.3|) 



2aX a(x + *) 

= — f + a?i 



I — a 2 1 — a 
we have (for convenience the index n is omitted) 



1 — or 



exp 



1.2(1 - a) J 



3 



(B.6) 



where the following notations are introduced 
i 

5 = J e- ta J Q [p s /z{l-z)y z , 
o 

X' =X+ (^-l)(x + t), 



aX 



<0, 



X 



> o. 



5(1 -a 2 ) 

Consider the integral 8- Using the definition of the Bessel function 

(-l) m {P/2) 2m {z - z 2 ) m 



(B.7) 

(B.8) 
(B.9) 



Jo [Pv z - z- 



E 

m=0 



m\ T(m + 1) 
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where T(x) is the Euler gamma-function, one represents the integral (1B.7P in 
the form 



3 = E 



m=0 



;-i) m (/3/2) 2 " 

ml T(m + 1) 



z - z 2 ) m dz 



(B.10) 



Calculating the latter integral [33fl, one obtains 



1 /2 



:-/3 2 /4i^ir 



ra=0 



'm+1/2 



M 

2 



(B.ll) 



where /^(x) is the modified Bessel function of the first kind. Furthermore, 
we use the equality jjyj 



t 



ml 



m=0 



Z 



-1/4 _ 



J 1/2 ( V^+2^ 



In our case 
M - [2*| 



z| - \2t\ > 0. 



(B.12) 



2a(l -a 2 

and instead of equation ( IB. Ill) we have 



a 2 7 2 + 2cry + a 2 ) > 



7T 



Z/|/i| 



- M /2j 



1/2 



where 



v = v\ 



S 2 



(5) = Jl-—(l- a 2 ) = - v / aV + 2a 7 + a 2 



a 



(B.13) 



(B.14) 



Substituting the expression f!B.13j) into equation (IB . 6|) and taking into ac- 
count that I\/<i(x) = (2/ixx) l l 2 sinh(x), we obtain 



Ji l) 



cxp 



X' 



2(1 - a) 



x 



x < exp 



ail - ^n)X 
'25 n (l-a 2 ) 



exp 



q(l + i/ ra )JT 
'25 n {l-a 2 ) 



, (B.15) 
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where v n = u(5 n ). At last, substituting X' from Eq. (1B.8[) into Eq. (1B.15P 
one has 



j(l) _ ^ e -7nX/2 



cxp 



2(1 - a 2 ) 



exp 



2(1 -a 2 



(B.16) 



and 



7, 



^ = a + 7 n ± a/ " 2 7 2 + 2 «7n + a 2 



r(2) 



TTie calculation of J„ 

Consider the integral JA . After substitution of the variable in equation 

dH3D 



X 



x — Xi 



1-a 2 

we have (the index n is omitted) 



(B.17) 



J(2) 



exp 



X' 



aX 



4a ^12(1 -a) 25(1 - a 2 



x 



J; 



.,f 1 Jl-a i If aX \2 
1-a 5 



01 2 V a 2 



1-a 2 ( aX \2 



(B.18) 



1-a 5 



To calculate the integral (IB. 18[) we consider the equality (1B.15P , having previ- 
ously made the substitution (IB.17P into Jn \ After reducing common factors, 
we have 



d^ 2S J 



1 /1-a 5 



a z 



aX \2 



1 — a 5 



sinh 



avX 



v L25(l-a 5 



(B.19) 
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Differentiating the latter equation with respect to X, one obtains 



X 
4a 



£/25_ 



J 2 



1 -« 2 



aX \2 



1-a 2 



^ 2 



a' 



aX \2 



1 — OZ' 



cosh 



cosh 



25(1 -a 2 ) 

One multiplies the latter equality by 

r X' aX ' 

~ 6XP 1.2(1 - a) ~ 25(1 -a 2 ). 

and using equations (IB. 18|) . (IB. 5|) and (IB.8|) . one has 

ln +) X 



avX 
25(1 - a 2 



J " = r l exp 27r^ +exp 27r^y} 



2 ieXP 



«7n(^ + t) + X 
2(1 - a) 



+ exp 



«7„(x — t) + X 
2(1 + a) 



(B.20) 



(B.21) 



Finally, substitute equations (IB.16j) and (IB.21[) into equation (IB.2|) and as a 
result we have 



j n (xt) = e- 7 " a ' /2 {^-)exp 



li +) X 
L2(l-a 2 



exp 



L2(l-« 2 



(B.22) 



where 



(±) _ An R 5n 



Appendix C. The determination of the parameters of the equa- 
tions (I33l)-(|36l) 

For the determination parameters entering (!33|) - (l36|) we use the condition 
continuity of the temperature across the interface, Ti = T$, and the condition 
of the heat balance ( IT2l) . For this purpose initially we write down the latent 
heat of solidification Q = kT^ with Tj = Ti\x=o = Tg\x=o, in dimensionless 
form as 

Q = Q/Q m = l + bf s \x= , (C.l) 
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where b = Q m /pc p T m . 

Now we equate the temperatures Tg and T L at the interface and substitute 
the fluxes fl34|) . fl36|) at X = into the condition of the heat balance f|T2|) . 
then taking into account the equalities f )26|) . one obtains 

ba ( s) + A /a = -1 (C.2) 
4 S) + A /a = A (C.3) 

Af = £rS (n £ 1) (C.4) 

—^±L + £^L_ = („ ;> 1) (C.5) 

n + 1 n n - b 

From the last equation of ( 1261) and the equality ( 1211) it follows that 

A^ — A^+i = A n . (C.6) 
The substitution of this equality into ( 1C.5I) gives the recurrent relationship 

^ +1 "n(2n+l-6)^ ^ lj ' 
whence one obtains 

^' = - (3-t)(5- t )...(2n-l- t) ^ «"» 2) ' (G7) 

where the equality A^ = —A has been used (see the relationships fl26|) ). 
The remaining parameters Aq and Oq are found from the solution of the 
system (1C.2P and (IC.3P in the form 

a 1 + bA (S) 1 + A 

Finally, taking into account the equalities ( 1C.4|) . ( 1C.7|) and (lC.8p . the expres- 
sions for the temperature of both the liquid and solid phase can be presented 
in the dimensional coordinates (x, t) in the form of the equations (1371) and 
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